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Calibration  of  the  HI  STAR  Sensors 


1.  INTRODUCTION 

The  .survey  data  compiled  into  the  "AFGL,  Four  Color  Infrared  Sky  Survey" 
catalog1  were  calibrated  against  selected  stars  observed  during  the  course  of  the 
scan  and  augmented,  when  necessary,  with  relative  res  pons  ivitles  measured  in  the 
laboratory.  The  celestial  sources  used  for  calibration  were  chosen  from  a list, 
compiled  at  AFGL,  of  objects  with  infrared  observations  reported  in  the  astronom- 
ical literature  prior  to  mid  1974.  Subjective  judgement  was  used  in  selecting  the 
most  reliable  sources  and  measurements  from  the  list. 

A detailed  description  of  the  experimental  profile  and  Instrumentation  is  given 
by  Price  and  Walker1  for  the  nine  rocket  flights  on  which  survey  data  were  obtained. 
Briefly,  the  survey  used  small  (16.  5 cm  diameter)  doubly  folded  Gregorian  tele- 
scopes cooled  to  cryogenic  temperatures.  The  focal  plane  of  each  telescope  con- 
tained three  linear  staggered  arrays  of  eight  detectors  each  arranged  in  the  cross 
scan  direction  and  spectrally  band  limited  with  interference  filters.  The  sensor 
scan  resulted  in  sequential  observations  in  the  three  spectral  colors.  The  cross 


(Received  for  publication  3 July  1978) 

1.  Price,  S.  D. , and  Walker,  R.  G.  (1976)  The  AFGL  Four  Color  Sky  Survey; 

Catalog  of  Observations  at  4.  2.  11.0.  lg, ,8  and  27.  4 um.  A FdL-tft -76-0205. 
Environmental  Research  Papers,  No.  578. 
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scan  neld-of-view  of  adjacent  detectors  In  a single  color  band  overlapped  at  least 
optical  one  blur  diameter. 

The  telescopes  used  for  the  rtrst  seven  experiments  were  filtered  for  the  4. 

11  and  20  micron  spectral  regions  and  used  doped  Germanium  detectors.  Silicon 
detectors  were  substituted  for  the  last  two  experiments  and  a spectral  band  cen- 
tered at  27  microns  was  used  In  place  of  the  4 micron  color.  The  same  sensor 
system  was  flown  on  the  first  six  experiments  and  different  telescope  and/or  focal 

plane  was  employed  for  the  last  three. 

A linear  regression  with  fixed  unit  slope  of  the  observed  system  magnitudes 
against  the  reference  magnitudes  was  performed  for  each  of  the  24  detectors  on 
each  flight.  The  preliminary  calibrations  on  the  first  six  experiments,  which  used 
the  same  sensor,  indicated  that  data  from  the  separate  flights  could  safely  be  com- 
bined to  Improve  the  accuracy  of  the  calculations.  This  sensor  was  also  twice 
calibrated  in  the  MARK  7 chamber  at  Arnold  Engineering  Development  Center 
(AEDC).  The  first  test2  was  performed  before  the  survey  experiments  began  and 
the  second3  was  done  after  the  sixth  flight.  Inconslstances  In  the  data  from  the 
first  test  obviated  the  use  of  these  results.  The  data  from  the  second  test  indicates 
the  presence  of  systematic  errors  in  the  aperature  sizes  assumed  by  AEDC  Tor  the 
blackbody  integrating  cavity.  However,  the  excellent  agreement  between  the  sec- 
ond AEDC  calibration  and  that  obtained  Trom  the  stars  indicates  that  these  errors 
were  small. 

There  were  an  Insufficient  number  of  reference  stars  observed  to  reliably 
calibrate  each  of  the  detectors  for  each  of  the  remaining  three  Rights.  The  refer- 
ence list  for  these  Mights  was.  therefore,  augmented  with  the  results  reduced  Trom 
the  previous  experiments.  Even  this  was  inadequate  In  some  cases,  especially  Tor 
the  27  Mm  color.  For  these  channels,  the  relative  res  pons  ivities  of  the  detectors 
In  each  color  measured  In  the  laboratory  for  the  focal  plane  alon  <■  i to 

the  available  celestial  reference  values.  The  internal  consistency  4ur« 

was  found  to  be  good. 


Hickey,  R.  F.  (1971)  Performance  Evaluation  of  the  SAM9Q /Hughes  HI  STAR 
Sensor.  AEDC-TR-71-B3  (Classified  Report). 

Nutt.  K.W.  (1973)  Performance  Evaluation  of  the  HI  STAR  Sensor  Ualnga 

New  Integrating  Sphere  Target  Source.  AEDC-TH-73-142  (Classified  Heport). 
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output  of  a linear  system  due  to  a spectral  irrndiance  Hx  is  given 


SXHX  dX  (1) 

where  Sx  is  the  system  spectral  response  and  R its  responslvlty  in  terms  of  volt/ 
watt.  These  quantities  are  routinely  measured  in  the  laboratory  for  the  focal  plane 
assembly  alone.  Invariably,  the  laboratory  calibration  sources  are  chopped  black- 
bodies  of  "known"  temperature  which  flood  the  full  field  of  the  detectors.  Thus, 
not  only  do  opt  teal  and  electronic  transfer  functions  have  to  be  adopted  to  convert 
the  measured  focal  plane  response  to  a system  response,  but  also  an  equivalence 
must  be  assumed  between  the  chopping  and  scanning  frequencies  and  the  illumination 
geometry  of  the  laboratory  sources  and  those  of  scanned  point-like  celestial  sources. 
Sayre  et  al4  have  measured  large  variation  in  responslvlty  with  illumination  geom- 
etry for  the  types  of  detectors  used  for  the  AFGL  survey. 

The  stars  offer  several  advantages  as  infrared  calibration  sources.  The 
entire  system  is  calibrated  during  data  taking  with  the  actual  point  source  scan 
geometry  and  with  sources  which,  for  the  majority  of  cases,  have  spectral  charac- 
teristics similar  to  that  of  the  cataloged  object.  Variations  in  responsivity  over 
the  surface  of  a detector  are  averaged  out  as  the  stellar  transit  positions  are  pot 
fixed.  Further,  for  practical  considerations,  stars  provide  the  only  reliable  cali- 
bration for  spectral  bands  such  as  the  4 micron  band.  This  filter  has  a two  micron 
wide  spectral  leak  of  about  4 percent  peak  transmission  at  14  pm.  Below  250°K  the 
blackbody  energy  through  this  "red"  leak  is  at  least  as  great  as  the  flux  in  the  pass- 
band.  Thus,  small  errors  in  the  detailed  knowledge  of  the  transmission  in  this 

leak  produce  large  errors  in  the  correction  terms  required  at  these  temperatures. 

2 3 

Unfortunately,  the  AEDC  * calibrations  used  blackbodies  with  temperatures  on  the 
order  of  2f>0°K.  The  known,  and  measured,  higher  temperature  of  the  reference 
stars  reduce  the  amount  of  the  red  leak  correction  for  these  sources  to  less  than 
one  percent. 

On  the  other  hand,  few  stars  have  been  measured  with  sufficient  spectral  detail 
in  the  infrared  to  allow  an  exact  evaluation  of  the  Integral  in  Eq.  (1).  Just  as  for 
the  laboratory  sources,  some  assumptions  are  necessary  concerning  their  spectral 


4.  Sayre,  C.,  Arrington,  D.,  Eisenman,  W.,  and  Merrlam,  J.  (1976)  Charac- 
teristics of  Detectors  Having  Partially  Illuminated  Sensitive  Areas,  IRIS 
Meeting  on  Detectors. 
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distributions.  Also,  there  are  too  few  true  astronomical  photometric  standard 
stars  bright  enough  to  reliably  calibrate  each  of  the  detectors  for  each  experiment, 
so  measurements  on  stars  of  late  spectral  type,  that  is,  cool,  small  amplitude 
ind  Mira  type  variable  stars  had  to  be  included  in  the  list  of  reference  sources. 

The  variability  of  these  stars  is  much  smaller  in  the  infrared  than  the  visual,  often 
amounting  to  only  10  to  20  percent. 

The  spectral  characteristics  of  a system  may  be  expressed  by  the  system 
unique  parameters  of  effective  wavelength,  designated  X and  effective  bandwidth. 
These  quantities  are  defined  by 


(2a) 


and 


AXe 


00 


(2b) 


respectively.  Table  1 lists  these  quantities  for  the  instruments  used  for  the  AFGL 
Infrared  Sky  Survey. 


Table  1.  Effective  Wavelength  and  Bandwidths  for  the  Survey 
Sensors 


■ ■ — 

Hi  Star 

Hi  Star  South 

(pm) 

&Ae  (pm) 

Xg  (pm) 

(pm) 

4.  16 

1.  50 

11.  11 

5.  67 

11.00 

5. 14 

19.63 

5.99 

19.  80 

5.  59 

27.43 

3.44 
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Expanding  the  irradiance  in  Eq.  (1)  about  the  effective  wavelength,  we  get' 

'‘-V2 <*-»/ 


VR'1  » HXe+J^~T 
Ae  2 dX2  x 


x /n'v|  / 


dX  (3) 


go 

/n‘V7 


<x  . xe)n  sx  dX V j sx 


and  j " o from  the  definition  of  *e. 

Equation  (4)  may  be  written 

VR  l " 1!Xe  A\»  (1  + AH)  <4) 

where  AH  is  the  color  correction  term  which  involves  the  second  and  higher  order 
terms  in  the  expansion.  This  is  an  exact  expression  equating  the  integrated  in- 
band  irradiance,  and  hence  the  system  response,  to  the  monochromatic  irradiance 
at  the  effective  wavelength  of  the  spectral  band  of  the  system,  H^e,  times  the 
respective  bandwidth  through  a proportionality  constant  AH,  a color  correction  term. 

_4 

For  the  spectral  bands  used  on  the  AFGL  Survey,  A11  is  about  0.25  for  a X infra- 
red spectral  distribution  and  less  for  flatter  variations  with  wavelength. 

The  photometry  on  the  reference  sources  is  in  terms  of  magnitude,  a logarith- 
mic relationship  between  the  stellar  irradiance  and  that  from  a standard  source. 

The  adopted  standard  source  is  an  AOV  star  at  a distance  of  10  parsecs 
19 

(3.  1 X 10  cm).  The  infrared  energy  distribution  from  such  a star  can  be  accu- 
rately represented  by  a 10,  000°K  blackbody  radiator  which  subtends  1.  5697  X 10’*b 
steradians.  The  magnitude  of  a star  at  X^  is  defined  by 


5.  Stromgren,  B.  (1937)  Handbuch  der  Experimental  Phvsik,  ed.  Wein  and  Harms. 
Leipzig:  Akademlshe  Verlags  Gesellschaft  M.B.  H.,  26:321. 
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mg(Xe>  = -2.  5 log 


HXe(a) 

HXe(AOV) 


The  observed  system  magnitude  is  defined  as 


n 


mQ  - -2.  5 log 


00 

J lysis* dA 

o 

oc 

J Hx(AOV)  Sj^  dX 


(5) 


which  can  be  expressed  in  terms  of  the  responsivity  from  Eq.  (1)  as 
* -2.  5 log  VR*1  + C 

o 

Here  CQ  is  the  logarithm  of  the  integrated  zero  point  flux  of  the  system  and  is 
uniquely  determined  for  each  color. 

Employing  the  relationship  given  in  Eq.  (4),  Eq.  (5)  can  also  be  rewritten 


) HXp(s)  AXp  (1  + AH(s)) 
m0  = "2-5  l°8  j HXe(AOV)  Me  (1  + AH(AOV)) 


(6) 


m = -2.  5 log 
o 


Hx  (s) 
HXe(AOV) 


-2.  5 log 


1 + AH(s) 

1 + AH(AOV) 


m = m (X  ) + C. 
o s e 1 


Equating  (6)  and  (8)  and  rearranging  terms,  we  have 


2.  5 log  V = -mg<xe)  + 2.  5 log  R + CQ  - Cj  . 


(7) 


(8) 


(9) 


V is  the  measured  variable,  voltage,  mg(X^)  is  the  reference  magnitude  input,  R 
is  the  sought  for  system  responsivity,  CQ  is  a known  system  value  and  is  a con- 
stant which  contains  an  unknown  color  correction  term.  Noting  that 


10 


Co  - Cj  ~ 2.  5 log  j HXe(AOV)  Me  + AH(s) 


if  AH(s)  < 0.  3 . (10) 

Equation  (9)  may  be  rewritten  in  terms  of  the  monochromatic  zero  point  flux: 

2.  5 log  V - -mg(*e)  + 2.  5 log  R + 2.  5 log  j HXe(AOV)  j + AH(s)  . (11) 

A reference  star  with  a known  magnitude  at  can  be  used  to  determine  the 
system  responsivtty  to  within  a constant  AH  by  means  of  Eq.  (11).  If  the  stars 
used  in  the  calibration  have  similar  infrared  energy  distributions,  then  their  color 
corrections  terms  are  approximately  the  same  and  average  calibration  may  be  ob- 
tained to  within  a mean  color  correction  value  representative  of  the  class  of  objects 
used.  The  star  to  star  differences  in  color  correction  terms  then  Introduces 
scatter  into  the  calculations. 

Almost  all  the  reference  stars  are  cool  giant  or  supergiant  stars  with  photo- 

spheric  effective  temperatures  between  2000  and  3500°K.  Many  of  these  stars 

possess  circumstellar  emission  which  enhance  the  flux  in  the  ten  micron  region  by 

as  much  as  a factor  of  three.  For  "oxygen"  stars  the  emission  feature  is  due  to 

silicate  grains.  In  "carbon"  stars  the  enhancement  is  caused  by  silicon  carbide 

plus  an  underlying  continuum  with  characteristic  temperatures  on  the  order  of 

600°K.  Detailed  studies  by  Merrill  and  Forrest  et  al  on  a number  of  the  objects 

in  the  reference  list  show  that  the  infrared  energy  distributions  for  these  stars  are 

smooth  and  somewhat  similar,  and  that  over  the  spectral  bands  used  for  the  survey 

the  energy  distributions  are  either  proportional  to  X for  the  photospheric  radiators 

or  somewhat  flatter  if  circumstellar  emission  is  contained  in  the  band.  Thus,  for 

the  system  bands  AH  varies  from  zero  for  an  equal  energy  distribution  to  about 
-4 

0.25  for  X distribution.  This  implies  that  the  approximation  used  in  Eq.  (10)  is 

6.  Merrill,  K.M.  (1977)  Infrared  Observations  of  Late  Type  Stars,  Invited  Review 

Paper  at  I.A.V.  Symposium  No.  42. 

7.  Forrest,  W.  J.,  Gillett,  F.C.,  and  Stein,  W.A.  (1975)  Circumstellar  grain 

and  the  intrinsic  polarization  of  starlight,  Ap.  J.  195:423. 
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valid  for  the  reference  llat  atara  and  the  atar  to  atar  color  correction  dlfferencea 
will  be  amall  aecond  order  terma. 

It  now  remalna  to  tranafer  the  reference  magnitudes  from  the  astronomical 
literature  to  valuea  at  the  system  effective  wavelengths.  The  X^  (4.2  (im)  refer- 
ence magnitudes  were  obtained  by  either  interpolating  linearly  In  magnitude  between 
published  values  at  3.  5 pm  and  5 jim  or,  If  there  la  no  5 pm  reference  valuea  re- 
ported, extrapolating  the  2.  2 pm  and  3.  5 pm  reference  valuea.  This  procedure  is 
valid  If  the  stellar  spectral  distribution  obeys  a power  law  in  wavelength  through 
this  region,  an  applicable  assumption  for  the  reference  list  sources.  A amall 
ayatematic  error  amounting  to  a few  percent  for  some  of  the  atara  may  have  been 
Introduced  due  to  CO  absorption  at  4.  3 pm. 

The  various  ground  based  photometric  systems  which  reported  stellar  meas- 
urements in  the  8 to  24  micron  region  have  color  bands  with  effective  wavelengths 
at  or  very  close  to,  that  is,  within  10  percent,  the  X^  (11.0  pm)  and  X^  (19.8  pm) 
system  bands.  The  ground  based  measurements  were  adopted  directly  and  errors 
due  to  mismatch  in  the  effective  wavelength  and  bandwidths  between  the  ground 
based  and  our  space  borne  systems  were  assumed  to  introduce  random  scatter  into 

the  calibration.  The  X (27.4  pm)  calibration  references  will  be  discussed 
c 

separately. 

The  compiled  list  of  reference  sources  is  contained  in  Appendix  A along  with 
the  adopted  magnitudes  at  the  system  effective  wavelengths.  The  maximum  and 
minimum  brightness  which  have  been  reported  are  also  given.  In  most  cases,  the 
adopted  magnitude  is  a simple  mean  of  the  brightness  extrema. 

Although  known  infrared  variables  were  accepted  as  reference  sources  under 
penalty  of  increased  scatter,  stars  known  to  vary  by  more  than  0.  75  magnitudes 
(a  factor  of  two  in  brightness)  were  eliminated  from  the  calibration.  Hearn  size 
effects  were  at  least  partially  accounted  for  by  deleting  all  objects  known,  or 
measured  by  the  survey,  to  be  extended. 

The  reference  sources  used  in  the  calibrations  were  further  restricted  to  stars 
8 9 

in  the  IRC.  ■ This  was  done  to  increase  the  probability  that  the  reference  source 
was  Indeed  a star  with  no  unusual  Infrared  spectral  properties  or  beam  size  effects. 

A linear  least  squares  regression  with  fixed  unit  slope  of  the  observed  magni- 
tude, more  precisely  2.  5 times  the  logarithm  of  the  observed  voltage,  against  the 
reference  magnitudes  is  calculated.  The  solution  Intercept  at  >us(Xp)  - 0 leads 
directly  to  the  system  responsivlty  through  Eq.  (11). 


8.  Neugebauer,  O. , and  l.elghton,  R.B.  (1909)  Two  Micron  Sky  Survey,  A 

Preliminary  Cntalog,  NASA  SP-3047. 

9.  Neugebauer,  G.  (1971)  Unpublished  extension  to  the  IRC. 


'WWW 


1 


J 


12 


S.  RESULTS 

The  sensors  proved  to  be  very  stable.  The  response  of  the  detectors  to  low 
temperature, evaporated  thin  film.  Infrared  emitting  sources  pulsed  during  the 
experiment  remained  constant  to  within  10  percent  even  though  the  background 
varied  by  as  much  as  a factor  of  five.  Furthermore,  the  preliminary  calibrations 
for  the  first  six  experiments  which  employed  the  same  sensor  system  had  remark- 
able homogenlety.  This  agreement  is  shown  in  Figures  1 and  2.  The  reference 
star  magnitudes  are  plotted  against  observed  magnitudes  for  two  of  the  Xg  (4.  2 pm) 
channels  In  Figures  la  and  lb,  while  Figures  2a  and  2b  are  similar  plots  for  two 
of  the  X (11,0  pm)  channels.  The  different  symbols  on  these  plots  represent 
different  flights.  The  data  from  the  Tlrst  six  experiments  were  combined  In  order 
to  Improve  the  accuracy  of  the  calculations.  Note  also  that  the  reference  sources 
extend  over  a factor  of  100  (5  stellar  magnitudes  in  the  figure)  in  brightness. 

An  iterated  linear  least  squares  regression  with  fixed  unit  slope  was  calculated 
for  each  detector  of  the  logarithm  of  the  observed  voltage  against  the  reference 
magnitude  by  means  of  Eq.  (11).  The  Intercept  of  the  solution  leads  directly  to  the 
system  responsivlty  if  an  average  value  is  adopted  for  the  color  correction  term. 
The  solutions  for  each  channel  were  Iterated  by  rejecting  as  many  as  five  sources 
with  the  largest  deviations  exceeding  twice  the  standard  deviation  of  the  previous 
solution.  This  was  done  In  order  to  minimize  the  effects  of  any  large  deviations 
due  to  unknown  beam  size  effects,  a greater  variability  of  the  source  than  the 
reference  values  Indicate,  or  large  unknown  measurement  errors.  Also,  only 
measurements  on  reference  sources  which  exceed  a specified  slgnal-to-noise  ratio 
were  included  In  order  to  minimize  the  measurement  error. 

Table  2 lists  the  least  squares  calibration  parameters  based  on  stellar  refer- 
ences for  the  combined  six  experiment  data  compared  to  other  calibration  methods. 
The  first  column  lists  the  channel  number;  channels  1 through  8 are  filtered  for  X^ 

(11,0  pm),  the  X (4.2  pm)  channels  are  9 through  15  (channel  16  was  inoperative) 
c 

and  channels  17  through  24  are  Xp  (19.  8 pm)  colors.  The  number  of  reference  star 
observations  used  In  calibrating  each  channel  after  all  rejections  is  listed  in 
column  3.  As  expected,  the  number  of  reference  sources  available  for  calibration 
decreases  with  increasing  wavelength.  Indeed,  there  were  hardly  enough  sources 
in  the  X (19.8  pm)  band  to  permit  a meaningful  regression  calculation. 

As  a check  on  the  linearity  of  the  system,  a separate  regression  was  calculated 
with  the  slope  as  a free  parameter.  Column  3 gives  the  value  of  the  derived  slope 
and  column  4 lists  the  root  mean  square  deviation  of  the  fixed  unit  slope.  Over  half 
the  channels  have  slopes  within  one  rms  of  unity  while  all  but  one  channel  are  with- 
in two  standard  deviations.  This  argues  strongly  that  the  system  response  was 
indeed  linear. 
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Table  2.  Comparison  of  the  Calibration  Parameters  Obtained  With  Stellar 
References  and  Values  Derived  From  Independent  Sources 
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A^(H  pm) 

No.  or 
Star* 

Kw 

Slop* 

KUH 

Slop* 

rma 

o 

mas 

o(K  * S 

man 

14  ” 

H'* 
i a i:  i k ' 

Nnta*) 

11  *l 

(Wak*> 

It*1 

<nt«Hi*t 

1 

17 

. 397 

. 058 

.315 

2.  988 

. 144 

3.  274 

2.  890 

3.087 

2.907 

2 

IN 

, 959 

.087 

. 294 

5.  294 

. 23 

3.  949 

5.  530 

5.  280 

5.922 

3 

It 

t.ots 

.094 

. 372 

2.  M»!> 

. 218 

3.098 

2.  99.1 

2.  71 

2.  085 

4 

23 

1 . 034 

. 054 

.329 

1.978 

. 153 

1.517 

2.  540 
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37 

. 990 
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1.82 
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7 

22 
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8 
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10 
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41 
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. 08 

3.178 
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12 

71 

l . 005 

.042 
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.04 
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2. 4 73 

5.  544 

13 

47 
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.041 
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4.055 

.04 
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4.  244 
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14 

72 
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15 
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1.  898 
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17 
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2.  5.18 

2.  848 
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7 
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1. 783 

2.022 

19 

9 
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.118 
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1.  395 

1 . 10 

2.  857 

2.  471 

4. 181 

3.  28 

*0 

It 

. !M7 
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. 258 

2.881 

. 70 

2.447 
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2.31.1 

2.41 

21 

4 
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.415 

2.858 

t.  10 

2.  771 

2.  824 

2.  815 

2.  21 

22 

It 
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.357 

l . 702 

. 58 

2.  .182 

2.  27 

1.911 

2.  28 

2.1 

7 

. 800 

.105 

. 51 1 

2.  828 

l.  19 

2.422 

2.  849 

4.12 

2.  79 

24 

9 

. 881 

. 094 

. 274 

2.420 

. 41 

2.  054 

1. 887 

1.72 

2.47 

The  standard  deviations  of  the  observed  and  calculated  magnitudes  are  given 
in  column  five.  These  values  are  not  significantly  different  than  the  variations  in 
responslvity  over  the  surfaces  of  the  detectors  measured  by  AKDC.  ’ ‘ 

The  inverse  responslvity,  It'1  in  Kq.  (1),  obtained  with  the  stellar  calibration 
references  and  normalized  to  unit  power  of  ten  is  given  in  column  <i  of  Table  2 and 
the  error,  in  terms  of  magnitude,  in  deriving  these  values  are  listed  in  column  7. 
For  small  values  of  error,  for  example,  o(H_1)  v 0.3,  the  magnitude  error  is 
approximately  equal  to  the  fractional  error  in  It"1,  that  is,  the  inverse  responslvity 
is  R 1 (1  ± <t(R  M).  For  comparison  the  inverse  res  pons  Ivttles  obtained  from  AKDC 

3 

tests  are  listed  In  columns  H and  9.  The  values  In  column  8 were  obtained  by  a 
linear  least  square  fit  of  the  measured  voltage  to  a "known"  Irradiance  from  a point 
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source  trails itting  a detector  at  the  nominal  survey  rate.  The  source  always 

transited  the  detectors  at  the  same  specified  elevation  so  these  values  are  not 

3 10 

averaged  over  the  detector's  surface.  The  AEDC  data'  ’ indicate  that  there  may 

be  a small  systematic  error  oT  5 to  15  percent  in  the  adopted  values  of  the  blackbody 

aperuture  sizes.  The  stundurd  deviations  of  the  observations  to  the  linear  fit  for 

the  X (11  pm)  values  in  column  8 are  between  (i  and  3(>  percent  and  are  comparable 
® "1 

to  the  equivalent  values  of  o(R  ) based  on  the  stars.  The  X (4.  2 pn\)  AEDC  data 

C 

have  large  discrepancies  which  are  certainly  due  to  the  red  leak  in  this  filter. 

The  inverse  responsivities  in  column  9 were  obtained  at  AEDC  by  measuring 
the  mean  square  noise  voltage  as  a function  of  a controlled,  known  background. 

The  measurement  errors  are  bound  to  be  somewhat  high  as  the  measured  voltage 
changes  with  the  square  root  of  the  photon  flux,  making  the  measurement  insensitive 
to  changes  in  the  input  variable  and  restricting  the  measurements  to  only  three  or 
four  points  because  of  dynamic  range  considerations.  Despite  these  factors,  the 
two  calibrations  are  within  30  percent  of  each  other  for  all  but  three  channels. 

Again,  the  red  leak  in  the  X^  (4.2  pm)  channel  precluded  meaningful  measurements 
by  this  technique. 

The  last  two  columns  contain  relative  inverse  responsivities,  scaled  to  the 

values  obtained  from  the  stars,  derived  from  transient  phenomena  which  fill  the 

field-of-view  of  each  detector  and  were  observed  on  each  flight.  These  values 

were  included  in  order  to  evaluate  possible  independent  measurement  methods 

which  could  be  used  to  derive  the  relative  detector  responses  and  to  check  possible 

geometry  effects.  Column  10  lists  the  relative  Inverse  responsivities  derived  from 

observations  made  on  the  high  altitude  re-entry  wake  of  the  rocket  sustalner.  As 

this  phenomenon  was  of  significant  angular  extent,  the  desired  information  was 

derived  by  extending  the  system  response  toward  dc  by  applying  the  inverse  bilinear 

z transform  or  the  electronic  high  pass  filter  to  the  data  stream.  The  resulting 

spatial  energy  distributions  were  well  represented  by  Gaussian  profiles.  The 

internal  consistency  of  the  resulting  amplitudes  were  between  5 and  13  percent  of 

the  scaled  mean.  The  amplitude,  flux  and  width  of  this  phenomena  are  observed  to 

change  from  scan  to  scan  over  a period  of  nine  seconds,  11  and  a time  dependence 

is  expected  during  a single  observation  in  the  sense  that  the  bottom  detectors  of  the 

array  observe  excited  air  from  a shock  which  was  formed  a fraction  of  a second 

more  recently  than  the  top  ones.  This  effect  cannot  be  large  as  the  measured  half 

widths  are  well  within  the  measurement  error  for  all  channels  on  a single  scan  and 

no  systematic  variation  is  evident.  The  X (11  pm)  and  X (19.8  pm)  inverse 

£ * 

10.  Hiatt,  J.  (1975)  A It  PA /Perkin -Elmer  MARK  VU  Sensor  Test,  AEDC-TR-75-50 

(Classified  Report). 

11.  Murdock,  T.  L. . and  Walker,  R.G.  ( 1975)  Infrared  Signatures  of  High  Altitude 

Reentry  Wakes.  AFCRI.-TR-75-0083  (Classified  Report). 
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responsivities  derived  from  the  wake  scale  to  within  a maximum  difference  of 
20  percent  or  the  stellar  calibrations.  The  Ae  (4.  2 pni)  discrepancies  are  probably 
due  to  the  red  leak  and  the  basic  nature  of  the  phenomenon.  The  physics  of  these 
events  have  been  discussed  in  more  detail  by  Murdock  and  Walker.  ^ 

The  scaled  inverse  responsivities  derived  from  the  internal  calibrators  are 
given  in  column  11.  Large  scale  geometry  effects  were  anticipated  due  to  the 
location  of  these  sources  on  one  side  of  the  focal  plane  housing.  These  measure- 
ments were  considered  on  the  assumption  that  the  internal  scattering  of  the  radia- 
tion within  the  focal  plane  housing  would  average  these  effects  out.  The  internal 
consistencies  of  these  pulses  were  also  good,  with  typical  deviations  of  20  percent 
of  the  scaled  mean  value.  These  relative  values  are,  for  the  most  part,  within 
30  percent  of  those  derived  from  the  stars.  Creator  discrepancies  exist  for  the 
A (4.  2 *im)  again  due,  in  part,  to  the  red  leak. 

Once  the  calibrated  inverse  responsivities  were  derived,  an  in  band  irradiance 

or  a magnitude  was  calculated  for  each  observation  by  means  of  Eq.  (1)  or  (fi) 

respectively.  The  A^  (4.  2 pm)  values  were  adjusted  for  the  red  leak  by  applying  a 

correction  term  derived  by  assuming  that  the  source  had  a blackbody  distribution 

in  thts  spectral  region  with  a color  temperature  defined  by  the  A^  (11.0  pm)  and 

A (19.8  pm)  observations.  The  correction  term  is  accurate  for  temperatures 
e 

greater  than  400  K.  Below  this  the  A (4.  2 pm)  measurements  are  very  uncertain. 

The  calibration  of  the  data  from  the  last  three  flights  were  more  difficult  as 
each  of  these  experiments  used  a different  focal  plane  assembly  and  thus  had  to  be 
calibrated  individually.  These  experiments  detected  an  insufficient  number  of  ref- 
erence stars  on  each  detector  to  provide  an  adequate  regression  solution,  so  the 
reference  list  was  augmented  for  each  subsequent  experiment  with  survey  measure- 
ments in  the  mutually  overlapping  area  fbom  previous  experiments.  The  calibra- 
tion errors  for  these  experiments  were  higher  due  to  the  larger  uncertainties 
inherent  in  the  previous  survey  measurements  and  the  fewer  total  number  of  sources 
to  calibrate  against. 

The  small  difference  in  wavelengths  between  the  11  pm  and  20  pm  shown  in 
Table  1 were  ignored  in  calibrating  the  southern  hemisphere  data.  The  A^  (ll.  l pm) 
and  A (19,  fi  wm)  colors  for  the  eighth  experiment  had  a sufficient  number  of  ref- 
erence sources  to  calibrate  each  detector  individually  by  the  previously  described 
technique. 

The  Inverse  responsivities  for  the  A^  (27.4)  channels  were  derived  by  scaling 

the  voltage  responsivities  measured  at  the  Naval  Electronics  Laboratory  Center 
1 2 

(NELC)  for  the  focal  plane  assembly  to  the  fluxes  inferred  from  published  data 


12.  Arrington,  D.C.,  and  Clsenman,  W.  L.  (1973)  Test  Data  for  a 24-Element 
Array.  2R00-15,  NELC.  
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on  GL2688,  GI.2495  = IRC  + 30407,  GL4114  n Car  and  GL2390  IRC  + 10420  and 
a blackbody  extrapolation  to  the  asteroid  Ceres  based  upon  the  color  temperature 
defined  by  the  X (11.  l^m)  and  X^  (19.  (i  pm)  observations.  The  estimated  calibra- 
tion errors  are  between  30  and  50  percent. 

The  scaling  factor  obtained  by  this  method  agrees  well  with  those  derived  from 
the  Xe  (11.  1 pm)  and  X^  (19.  6 pm)  individual  channels  and  is  essentially  the  total 

system  gain  for  each  color.  Unfortunately,  these  system  gains  were  4 to  4.  5 times 

13 

less  than  those  calculated  from  the  known  optical  gain  of  the  telescope  and  the 
electronic  amplification  as  measured  in  the  laboratory.  If  the  amplifier  gains  are 
applied  to  the  detector  noise  values  measured  at  NELC,  the  result  is  within  50  per- 
cent of  that  observed  during  the  flight. 

The  measured  lower  responsivity  of  the  detectors  was  originally  attributed  to 

4 

the  illumination  geometry  effects  observed  by  Sayre  et  al  for  parallel  biased 

12  14 

detectors  such  as  the  ones  used  for  the  survey.  However,  Arrington  et  al  ’ 
noted  that  the  detector  responsivities  decreased  with  increasing  background  flux. 

The  NELC  data^'  **  indicated  that,  if  account  is  made  for  the  "blue"  spectral  leak 
in  the  X (27.4  pm),  the  character  of  the  responsivity  decrease  in  each  of  the 
colors  could  be  explained  by  a dilute  300°K  background.  This  proved  to  be  the  case 
as  an  increase  of  the  proper  amount  was  observed  in  the  detector  response  to  the 
internal  calibrators  when  the  focal  plane  was  optically  sealed  off  from  the  rest  of 
the  sensor.  This  decrease  in  sensor  sensivity  unfortunately  about  compensated 
for  the  anticipated  increase  due  to  the  sensor  modifications  performed  for  the 
southern  hemisphere  experiments. 

After  the  calibrations  were  applied  to  the  individual  survey  experiments,  j 

multiple  observations  on  a single  object  obtained  from  the  different  flights  were  j 

combined  in  a simple  mean.  Comparisons  between  the  resulting  final  survey  mag-  1 

nitudes  and  reference  source  magnitudes  for  all  objects  common  to  both  lists  are 

shown  in  Figures  3,  4,  and  5.  In  Figure  3 the  difference  between  the  X^  (4.2  pm) 

measured  and  reference  magnitudes  are  plotted  against  the  measured  magnitudes. 

The  vertical  lines  in  this  plot  are  the  reported  extreme  on  the  reference  sources. 

The  same  plot  for  the  X (11.  0 pm)  color  is  shown  in  Figure  4.  Here  the  crosses 
e 15 

are  magnitudes  taken  from  the  work  of  Hall  who  adopted  a different  set  of  refer- 
ence values  for  many  of  the  sources  in  common  with  the  AFGL  list.  Finally,  the  ! 

same  data  for  the  X (19,  8 pm)  is  shown  in  Figure  5.  | 

- — 

13.  HI  STAR  II  PROGRAM  (1974)  Final  Report,  Hughes  Aircraft  Company,  Report 

No.  p74-288  (Classified  Report). 

14.  Arrington,  D.  C.,  and  Eisenman,  W.  L.  (1974)  Test  Data  for  a 24-Element 

Array,  2610-17,  Naval  Electronics  Research  l.abortilorlfts.  j 

15.  Hall,  R.T.  (1974)  A Catalog  of  10-pm  Celestial  Objects.  SAMSO-TR-74-212. 


MEASURED  MAGNITUDE 


Figure  3.  The  Measured  Minus  Reference  Star  Magnitudes  Plotted  Against  the 
Measured  Magnitude  at  Xe  (4.2  urn)  for  all  the  Reference  List  Sources  Observed 
During  the  AFGL  Survey.  The  vertical  lines  represent  the  range  of  values  re- 
ported in  the  astronomical  literature  while  the  symbol  bisected  by  the  line  is  the 
adopted  mean  reference  value  for  this  source 


The  AFGL  survey  is  estimated  to  be  statistically  complete  down  to  magnitudes 
of  m(4.2  (im)  i +1.3,  m(11.0jim)~  -1.  1 and  m(19.  8 |im)  ~ -3.1.  Below  these 
levels  the  measurement  error  is  a significant  factor  in  calculating  the  observed 
magnitudes.  This  is  evident  in  Figures  3,  4,  and  5 where  the  difference  between 
the  measured  and  reference  magnitudes  become  rather  large  for  measured  magni- 
tudes fainter  than  the  aforementioned  limits.  Also  it  may  be  noted  that  at  the 
fainter  magnitudes  there  is  a tendency  for  the  measured  quantities  to  be  brighter 
than  the  reference  values.  This  can  be  explained  by  the  fact  that  at  these  levels  a 
positive  measurement  error  added  to  the  true  observed  voltage  of  the  source  will 
yield  a brighter  than  correct  measurement  while  a negative  error  will  cause  the 
observation  to  fall  below  the  detection  threshold  and,  therefore,  eliminate  it  from 


MEASURED  MAGNITUDE 

Figure  4.  Comparison  of  the  Xe  (U.  0 pm)  Measured  Magnitudes  and  Their  Re- 
spective Measured  Minus  Reference  Magnitudes  for  all  Reference  List  Sources 
Observed  During  the  Survey.  The  points  dotted  with  x's  are  reference  magni- 
tudes taken  from  the  compilation  of  Hall. 15  Note  that  for  most  of  the  variable 
stars  the  value  adopted  by  Hall15  is  different  than  that  from  the  AFGL,  compila- 
tion. The  other  symbols  have  the  same  meaning  as  Figure  3 


consideration.  The  threshold  for  a measurement  to  be  included  in  the  calibration 
was  set  to  be  equivalent  to  the  limit  of  statistical  completeness  in  each  color,  so 
this  last  effect  should  have  been  minimized  in  the  calibration. 

An  indication  of  the  accuracy  of  the  calibration  may  be  obtained  from  Table  3 
which  compares  the  observed  magnitudes  to  the  adopted  standard  magnitudes  for 
28  bright  stars  which  are  well  measured  in  the  infrared  and  are  either  known  to  be 
non-variable  or  have  been  used  as  infrared  standard  stars  for  various  ground  based 
photometric  systems.  As  can  be  seen  the  two  magnitudes  agree  within  a few  tenths 
of  each  other.  For  the  stars  in  Table  3,  the  observed  values  in  each  color  are,  on 
the  average,  a tenth  of  a magnitude  brighter  than  the  standard  magnitudes  with 
root  mean  squared  deviations  of  0.2  magnitude  about  this  mean  difference. 


MEAS-REF  MAGNITUDE 


Figure  5.  Comparison  Between  the  Measured  and  Measured  Minus  Reference 
Magnitudes  at  (19.8  ^m)  for  all  Sources  in  Common  With  the  AFGL  Catalog1 
and  the  Reference  List.  Note  the  paucity  of  objects.  The  symbols  have  the 
same  meaning  as  in  Figure  3 
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Table  3.  Observed  and  Reference  Magnitudes  for  Stars  Commonly  Used  as 
Infrared  Standards 
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4.  CONCLUSIONS 


Several  significant  observations  may  be  made  from  the  results  of  the  calibra- 
tion procedures  used  at  AFGl.  for  its  space  borne  cryogenically  cooled  sensors: 

(1)  The  linear  slope  over  two  orders  in  magnitude  for  each  detector  obtained 
from  the  least  squares  regression  of  the  system  magnitudes  against  the  reference 
values  for  the  combined  data  from  the  first  six  flights,  attest  not  only  to  the 
linearity  and  stability  of  the  system  but  demonstrate  that  stellar  references  are 
good  calibration  sources  provided  enough  of  them  are  observed. 

(2)  Primary  standards,  that  is  sources  with  precisely  known  spectral  energy 
distributions,  are  not  required.  Spectral  uncertainties  and  amplitude  errors  for 
small  amplitude  variable  stars  add  in  a random  fashion  such  that  the  average  cali- 
bration is  correct.  If  enough  stellar  references  are  observed  over  an  adequate 
dynamic  range,  the  error  in  calculating  the  system  responsivity  is  smaller  than 
the  error  in  the  regression  fit.  Further,  the  standard  deviation  of  the  observed 
and  reference  magnitude  is  of  the  same  size  as  the  measured  variation  in  respon- 
sivity over  the  surface  of  a detector. 

(3)  Rather  surprisingly,  a good  calibration  was  obtained  from  a least  squares 
fit  of  the  square  of  the  observed  noise  voltage  to  a known  background  irradiance. 
Although,  this  procedure  suffers  from  lack  of  sensitivity  and,  for  practical  reasons 
restricts  the  dynamic  range  of  the  independent  variable,  it  may  have  value  in  cali- 
brating some  of  the  spectral  colors  proposed  for  future  systems  where  few,  if  any, 
reliable  laboratory  or  celestial  sources  are  known  with  any  accuracy. 

In  conclusion,  it  was  found  that  the  stars  provided  excellent  reference  sources 
to  calibrate  space  borne  sensors  if  a sufficient  number  of  them  are  observed  to 
permit  a meaningful  regression  solution. 
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Appendix  A 


The  fable  in  this  appendix  list  the  IRC  source  which  comprised  the  reference 
star  list  used  to  calibrate  the  AFGL  sensors.  The  first  column  lists  the  IRC  num- 
ber of  the  star  while  the  second  gives  the  star  name,  variable  star  or  other  desig- 
nation. The  maximum,  mean  and  minimum  magnitudes  at  (4,  2 pm)  adopted  for 
this  source  is  given  in  columns  3,  4 and  5,  respectively.  The  respective  A# 

(11.0  pm)  values  are  listed  in  columns  fi,  7 and  (i  while  columns  P,  10  and  11  give 
the  A^  (19.8  pm)  magnitudes. 

Note  that  in  the  name  column,  MIHJ  and  Nttl'  were  used  to  designate  p and  v 
in  order  to  avoid  possible  confusion  with  variable  star  designations  of  Ml'  and  NU. 
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